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Abstract 

A  two-dimensional  isothermal  model  of  an  SOFC  operating  with  a  syngas  fuel  based  on  button  cell  geometry  was  developed.  Momentum,  mass 
and  charge  transport  coupled  with  electrochemical  and  chemical  reactions  were  considered.  The  developed  model  has  been  validated  by  showing 
good  agreement  with  experimental  data  for  various  syngas  compositions  obtained  at  an  operating  temperature  of  900  °C.  However,  the  model 
predictions  deviate  somewhat  from  the  experimental  data  at  800  °C.  Visual  and  SEM  investigation  of  the  anode  surface  after  completion  of  each 
experiment  showed  that  excessive  amounts  of  carbon  form  on  the  anode  surface  and  that  the  Ni  catalyst  delaminates  from  the  anode  during  cell 
operation  at  800  °C.  In  contrast,  no  carbon  formation  was  observed  and  only  a  small  amount  of  Ni  delaminated  from  the  anode  during  operation 
at  900  °C.  The  model  was  further  applied  to  predict  carbon  formation  observed  experimentally,  and  to  study  the  impact  of  operating  conditions  on 
carbon  formation  and  determine  a  suitable  range  of  conditions  that  suppress  its  onset. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

With  the  rising  concern  about  greenhouse  gas  emissions, 
many  countries  are  pursuing  efforts  to  develop  more  sustainable 
energy  systems  to  replace  conventional  combustion  heat 
engines.  Solid  oxide  fuel  cell  (SOFC)  power  generation  shows 
great  promise  to  serve  as  an  alternative  in  the  near  future. 
Since  SOFC  operation  relies  on  electrochemical  reactions,  its 
electrical  efficiency  is  not  limited  by  temperature  as  in  the 
case  of  conventional  heat  engines.  Therefore,  the  electrical 
efficiency  obtained  from  a  SOFC  is  typically  greater  and 
more  stable  than  that  obtained  from  conventional  heat  engines 
[1,2].  Furthermore,  additional  efficiency  can  be  gained  by 
incorporating  a  steam/gas  turbine  cycle  to  recover  heat  from 
the  hot  gas  exhausted  from  SOFCs  because  they  operate  at  high 
temperatures  between  800  and  1000  °C  [1].  Fuel  flexibility  with 
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SOFCs  also  presents  advantages  over  other  types  of  fuel  cells 
such  as  polymer  electrolyte  membrane  fuel  cells  (PEMFCs). 
Not  only  is  the  carbon  monoxide  present  in  the  fuel  gas  not 
harmful  to  SOFCs,  but  it  can  also  act  as  a  fuel.  Therefore,  a 
variety  of  hydrocarbon-based  fuels  or  their  syngas  derivatives, 
such  as  natural  gas,  biomass  and  coal,  can  efficiently  be  used  in 
SOFCs. 

Biomass  and  coal  are  gasified  to  form  syngas,  essentially  con¬ 
taining  a  mixture  of  H2,  CO,  CO2,  H2O,  N2  and  eventually  H2S. 
This  syngas  can  then  be  used  as  a  fuel  in  SOFC.  However,  H2S 
must  be  scrubbed  before  the  syngas  is  fed  to  the  SOFC  to  prevent 
rapid  degradation  of  the  SOFC  performance  [1].  Depending  on 
the  type  of  coal  and  biomass,  the  composition  of  syngas  varies 
widely  and  affects  the  SOFC  performance  [3-8].  Also,  carbon 
formation  is  a  problem  that  can  limit  SOFC  performance  when 
used  with  syngas.  This  is  due  to  the  Boudouard  reaction  which 
converts  CO  to  carbon  on  the  anode  surface: 

2CO  ->  C  +  C02  (1) 
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Nomenclature 

c 

concentration  (mol  m-3) 

Dij 

binary  diffusivity  of  species  i  in  species  j  (m2  s-1) 

Di-M 

binary  diffusivity  of  species  i  in  multi-component 
gaseous  mixture  (m2  s-1) 

£° 

ideal  half-cell  cell  potential  (V) 

F 

Faraday’s  constant  (96487  A  s  mol-1) 

4vg 

average  current  density  (Am-2) 

J 

current  density  (A  m-2) 

4vgsr 

reaction  rate  constant  for  the  water  gas  shift  reac¬ 
tion 

Keq 

equilibrium  constant 

M 

molecular  weight  (g  mol-1) 

p 

pressure  (Pa) 

r 

radial  coordinate  (m) 

R 

gas  constant  (8.31441  Jmol-1  K-1) 

^contact 

contact  resistance  (£2  m-2) 

^WGSR 

flux  caused  by  WGSR  (molm-2  s-1) 

interface  conductivity  for  electrochemical  reac¬ 
tion  (£2-1  m-2) 

T 

operating  temperature  (K) 

U ,  V 

velocity  components  in  r  and  z  directions,  respec¬ 
tively  (ms-1) 

^delivered 

cell  potential  delivered  to  the  load  (V) 

Z 

axial  coordinate  (m) 

Greek  letters 

a  c 

carbon  activity  for  the  Boudouard  reaction 
defined  in  Eq.  (23) 

0 

overpotential  (V) 

M 

viscosity  (kgm-1  s-1) 

P 

gas  density  (kgm-3) 

ionic  conductivity  (£2-1  m-1) 

4> 

potential  (V) 

0an>  0cat 

half-cell  potential  at  the  anode  and  the  cathode, 
respectively  (V) 

Subscripts 

an 

anode 

B 

Boudouard  reaction 

cat 

cathode 

e 

electrolyte 

f 

fuel  channel 

i 

fuel  gas  species  (H2,  H2O,  CO,  CO2  and  N2) 

WGSR 

water-gas  shift  reaction 

Therefore,  various  operating  strategies  to  suppress  carbon  for¬ 
mation  must  be  devised.  Faced  with  these  potential  problems, 
it  is  important  to  develop  a  model  to  predict  SOFC  perfor¬ 
mance  for  various  operating  conditions,  in  particular,  syn¬ 
gas  feed  composition.  Such  a  model  could  also  be  used  to 
study  the  influence  of  these  operating  conditions  on  carbon 
formation. 


In  this  paper,  a  mechanistic  SOFC  transport  model  coupled 
with  electrochemical  (H2  oxidation,  CO  oxidation  and  O2  reduc¬ 
tion)  and  chemical  reaction  (water-gas  shift)  kinetics  for  syngas 
operation  (mixtures  of  H2,  H2O,  CO,  CO2  and  N2)  based  on  a 
button  cell  geometry  is  developed.  A  unique  advantage  of  this 
model  is  that  it  accounts  for  electrochemical  H2  and  CO  oxi¬ 
dation.  The  developed  model  is  validated  using  experimental 
cell  performance  data  obtained  for  various  syngas  composi¬ 
tions  at  800  and  900  °C.  Visual  and  SEM  investigation  of  the 
anode  surface  after  cell  operation  at  each  operating  tempera¬ 
ture  is  performed.  The  model  is  also  applied  to  predict  carbon 
formation  observed  during  experiments  and  then  employed  to 
study  the  influence  of  operating  conditions  on  carbon  deposi¬ 
tion. 

2.  Experimental 

The  electrochemical  cell  used  in  this  study  was  an  Electrolyte 
Supported  Cell  (ESC)  supplied  by  Innovative  Dutch  Electro 
Ceramics  InDEC  B .  V.  Co.  Ltd.  The  ESC  consists  of  a  TZ3  Y  elec¬ 
trolyte  (3  mol%  yttria- stabilized  zirconia)  that  is  approximately 
150  jmm  thick  and  fabricated  by  tape  casting.  A  60  p,m  thick 
Ni-Ce02-YSZ  anode  and  a  50  pan  thick  Lao.7Sro.2Mn03_<5 
(LSM)  cathode  were  screen  printed  on  each  side  of  the  elec¬ 
trolyte.  The  cross-sectional  areas  of  the  electrolyte  and  each 
electrode  are  25  cm2  (5  cm  x  5  cm)  and  16  cm2  (4  cm  x  4  cm), 
respectively. 

A  schematic  of  the  test  rig  is  illustrated  in  Fig.  1.  The  anode 
part  of  the  ESC  was  placed  on  a  Ni  mesh  current  collector  spot- 
welded  to  two  Ni  wires  acting  as  current  and  voltage  probes. 
Both  the  Ni  mesh  and  the  ESC  were  fixed  by  a  high  temper¬ 
ature  glass  seal  to  the  top  of  a  quartz  tube  which  was  in  turn 
attached  to  a  stainless  steel  sample  holder  by  a  ceramic  seal. 
The  diameter  of  the  quartz  tube  was  45  mm  i.d.  with  3  mm  thick¬ 
ness.  Two  smaller  quartz  tubes  (6  mm  diameter)  with  different 
lengths  were  fixed  to  the  sample  holder  using  high  tempera¬ 
ture  ceramic  seals.  The  long  tube  was  used  for  gas  inlet  while 
the  shorter  one  served  as  the  tube  for  gas  outlet.  The  gas  feed 
tube  was  located  4-5  mm  away  from  the  cell  to  ensure  that  the 
reactant  gas  was  uniformly  distributed  over  the  entire  area  of 
the  cell.  On  the  cathode  side,  another  quartz  tube  was  placed 
on  the  cell  to  press  the  glass  seal  when  it  melted.  This  ensured 
that  the  glass  seal  fused  to  the  cell  and  no  leak  occurred.  The 
cathode  part  was  directly  exposed  to  ambient  air  so  that  air 
flowed  to  the  cell  via  natural  convection.  A  platinum  mesh 
(Alfar  Aesar  Co.  Ltd.)  was  used  as  the  current  collector  at  the 
cathode.  A  weight  of  500  g  was  placed  on  top  of  the  platinum 
mesh  to  ensure  that  it  was  in  good  contact  with  the  cathode. 
This  also  helped  reduce  contact  resistance  between  the  plat¬ 
inum  mesh  and  the  ESC.  A  highly  porous  alumina  layer  was 
placed  between  the  weight  and  the  platinum  mesh  to  allow  the 
oxidant  (oxygen)  to  reach  the  cathode  reaction  site  underneath 
the  weight  uniformly.  The  cell  test  rig  was  contained  inside  a 
temperature-controlled  furnace  (Carbolite  Co.  Ltd.)  to  ensure 
that  it  could  be  operated  under  isothermal  conditions.  Two  ther¬ 
mocouples  were  located  at  the  vicinity  of  the  anode  and  the 
cathode. 
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Sling  for  weight 


Fig.  1.  Schematic  diagram  of  test  rig  set  up. 


During  the  operation  of  the  experimental  set-up,  H2,  CO,  CO2 
and  N2  were  mixed  together  using  a  set  of  mass  flow  controllers 
(Omega  Co.  Ltd.).  The  gas  mixture  was  then  passed  through  a 
flame  arrester  and  a  hydrator  prior  to  entering  the  anode  side  of 
the  electrochemical  cell.  A  hydrator  placed  on  a  temperature- 
controlled  bath  was  used  to  humidify  the  inlet  gas  stream.  The 
extent  of  gas  humidification  was  controlled  by  the  bath  tem¬ 
perature.  The  exit  gas  from  the  anode  side  was  vented  to  the 
atmosphere  through  a  fume  hood.  The  inlet  and  exit  lines  were 
connected  to  a  gas  chromatograph  (GC)  to  measure  gas  compo¬ 
sitions. 

The  cell  performance  was  characterized  by  measuring  the 
cell  potential  as  a  function  of  the  current  density  for  various  gas 
compositions,  gas  flow  rates  and  operating  temperatures.  The 
cell  potential  and  current  were  recorded  only  after  the  response 
had  maintained  a  steady  state  for  at  least  3  min  after  the  impo¬ 
sition  of  each  load  resistance  setting.  The  load  resistance  was 
varied  from  high  to  low  values  corresponding  to  a  range  from  low 
to  high  currents  drawn  from  the  cell.  In  these  experiments,  the 
influences  of  N2  and  CO2  dilution,  H2/CO  composition  ratio 
and  simulated  syngas  composition  (mixture  of  H2,  CO,  CO2 
and  N2)  on  cell  performance  were  investigated.  The  operat¬ 
ing  temperatures  for  these  experiments  were  800  and  900  °C. 
After  each  experiment,  the  anode  surface  was  inspected  visu¬ 
ally  and  by  scanning  electron  microscope  (SEM)  to  observe  any 


change  in  the  anode  microstructure  and  any  eventual  carbon  for¬ 
mation. 

3.  Mathematical  modeling 

3.1.  Model  geometry  and  assumptions 

The  geometry  of  the  test  rig  is  based  on  a  button  cell  design 
including  an  anode  feed  tube,  air  channel,  fuel  channel  and  elec¬ 
trochemical  cell.  Due  to  symmetry  conditions,  only  the  portion 
of  the  rig  shown  in  Fig.  2  is  considered  in  the  model  using  two- 
dimensional  cylindrical  coordinates  (radial  and  vertical  direc¬ 
tions  r  and  z).  Fuel  gas  is  supplied  through  the  anode  feed  tube. 
Once  it  reaches  the  anode,  it  flows  radially  outward  along  the 
porous  anode  layer  where  it  participates  in  the  electrochemical 
and  the  water-gas  shift  reaction  (WGSR)  reactions.  Gas  is  dis¬ 
charged  from  the  cell  by  flowing  from  the  anode  layer  into  the 
gas  outlet  channel.  Since  the  cathode  is  directly  exposed  to  the 
ambient  air  within  the  furnace,  no  feed  tube  is  required. 

Assumptions  used  in  this  model  development  are  as  follows: 

(a)  Temperature  is  uniform  over  the  entire  model  geometry. 
This  assumption  is  reasonable  because  the  test  rig  is  small 
and  is  placed  in  a  controlled-temperature  furnace.  Although 
heat  is  generated  by  the  electrochemical  reactions,  tem¬ 
perature  only  slightly  increases  since  the  electrode  area  is 
relatively  small  and  so  heat  generation  is  not  significant. 

(b)  Gas  density  is  constant  because  the  temperature  is  uniform 
and  the  generation  and  consumption  rates  of  gas  species  are 
low. 

(c)  Mass  transport  and  ohmic  resistances  within  the  anode  and 
the  cathode  are  neglected  since  their  thicknesses  are  very 
thin.  This  assumption  is  based  on  the  work  of  Chan  et  al.  [9] . 
Therefore,  electrode  subdomains  are  ignored  in  the  model. 
Only  the  ohmic  resistance  through  the  electrolyte  layer  is 
considered  in  the  model  because  the  ionic  conductivity  is 
much  less  than  the  electronic  conductivity  of  the  electrodes. 

(d)  The  electrochemical  reactions  and  the  WGSR  occur  exclu¬ 
sively  at  the  electrode/electrolyte  interfaces. 


r.  =  18  mm 

I) 

Fig.  2.  Model  geometry  of  the  test  rig  (not  to  scale). 
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Table  1 

List  of  governing  equationsa 


Subdomain  Transport  equation 


Fuel  channel 
Overall  mass 

Momentum  z  direction 

Momentum  r  direction 

Species 

Electrolyte 

Ion 

Air  channel 
Oxygen 


dvf  |  1  djruf)  _ 
dz  r  dr 
pd(vfVf)  1  pd(rufVf) 

dz  r  dr 

pd(vfUf)  1  pd(rufUf) 

dz  r  3r 

d(cjV f)  1  d(rciUf)  _  3 

dz  r  dr  dz 


dpf  3 

—  +  — 
dz  dz 

dpt  3 

—  H - 

3r  dz 


Di-M 


dcj 

dz 


Variable:  0e 
1  d(rd<t>e/dr )  ,  a20e 
r 


+ 


3r 

Variable:  co2 

a  /  dc0l 

-Id°2_N2  — 


dz2 


=  0 


dvf  \  13/  dvf 

Pf^~  +  -  V  (  Wv 
9z  /  r  dr  \  dr 

'  duf\  1  3  /  3wf' 

/v  J  +  r/if  — 

dz  /  r  dr  \  dr 

1  3  /  3c,- 

d —  TT  (  r Di-M  — 
r  dr  \  dr 


1  3  f  dc0l 

+  ~rJr  (rD02-N2lT 


=  0 


3  /  3uf 

+  —  ix  f  — 

3z  V  9z 

3  /  dvf 

+  9z  V  f_9r 


13/  3wf 

+  ^lrwaT. 

13/  3wf' 


2/ifUf 


(2) 

(3) 

(4) 

(5) 


(6) 


(7) 


a  Where  i  are  H2,  H2O,  CO  and  CO2,  Mf  and  Uf  are  the  fuel  velocities  in  r  and  z  directions,  respectively,  pf  the  pressure  field  within  fuel  channel,  ct  the  concentration 
of  ith  species,  Di-m  the  binary  diffusivity  of  ith  species  in  multi-component  mixture,  <pe  the  potential  within  the  electrolyte  subdomain  layer,  co 2  the  partial  pressure 
of  O2  and  Dq2  _n2  is  the  binary  diffusivity  of  O2-N2  system. 


3.2.  Governing  equations 


The  governing  equations  are  summarized  in  Table  1.  The 
Navier-Stokes  equations  [10]  and  mass  transport  equations  [11] 
within  the  fuel  channel  are  solved  simultaneously  to  determine 
the  velocity  of  the  fuel  stream  and  concentration  profile  of  each 

Table  2 


List  of  operating  conditions  and  model  input  parameters 


Item 

Value 

Operating  conditions 

Pressure  (Pa) 

1  atm  (101.35  kPa) 

Temperature 

800  and  900  °C 

Fuel  flow  rate 

100-180  ml  min-1 

Fuel  gas  composition 

Vary 

Cell  potential  (0cat  -  0an) 

0.5-1. 0  V 

Oxidant 

Air 

Model  input  parameters 

Conductivity  of  the  electrolyte  (cre, 

2.26  at  800  °C,  5. 12  at  900  °C 

f^m-1) 

Equilibrium  and  kinetic  constants 

Ke q,H2  (atm  0  5 ) 

1.78  x  109  at  800  °C, 

1.66  x  108  at  900  °C 

^eq,CO  (atm-0-5) 

1.63  x  109,  1.07  x  108 

^eq,WGSR 

1.04,  0.75 

£w gsr  (molm-2  atm-2  s-1) 

1.6  X  10s,  4.3  X  108 

Binary  diffusion  coefficients  (cm2  s  1 ) 

Du2-n2o 

7.68  at  800  °C,  8.94  at  900  °C 

Du2-  co 

6.45,  7.47 

Dh2-co2 

8.67,  10.0 

dh2-  n2 

6.36,  7.38 

Oh2o-co 

2.22,  2.58 

Dh2o-co2 

1.71,2.01 

Du2o-n2 

2.19,2.58 

Dco-co2 

1.41,  1.62 

Dqo-n2 

1.74,2.04 

Dc  o2-n2 

1.38,  1.62 

do2-n2 

1.83,2.07 

Fig.  3.  Experimental  (symbols)  and  predicted  (lines)  cell  performances  at  vari¬ 
ous  fuel  flow  rates  at  (a)  800  °C  and  (b)  900  °C  using  H2  as  fuel. 
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gas  species.  Fuel  gases  considered  in  the  model  are  H2,  CO, 
CO2,  H2O  and  N2.  Only  charge  transport  [12]  is  considered 
within  the  electrolyte  subdomain  layer  to  calculate  the  potential 
distribution  and  calculate  the  ohmic  overpotential.  Mass  trans¬ 
port  and  ohmic  resistances  within  the  porous  electrode  layers  are 
ignored.  Only  mass  transport  of  O2  is  considered  within  the  air 
channel.  Electrochemical  and  chemical  reactions  considered  in 
this  model  are  O2  reduction  at  the  cathode,  H2  and  CO  oxidation 
at  the  anode,  and  the  WGSR  at  the  anode/electrolyte  interface, 
which  are  described  in  detail  next. 

3.3.  Electrochemical  reactions 

Electrochemical  reactions  are  coupled  with  the  transport 
equations  to  determine  the  distributions  of  the  local  current 
density  and  electrode  overpotentials  inside  the  SOFC.  They  are 
specified  here  as  boundary  conditions  at  the  electrode/electrolyte 
interfaces.  In  reality,  electrochemical  reactions  in  an  SOFC  take 
place  at  the  triple  phase  boundaries  (TPB)  between  the  ionic 
conductor,  electronic  conductor  and  gas  in  the  vicinity  of  elec¬ 
trode/electrolyte  interface  where  O2-  and  electrons  meet  [13]. 
Important  aspects  and  the  rate  expressions  of  these  reactions  are 
summarized  in  this  section. 

3.3.1.  Oxygen  reduction  reaction 

The  rate  of  oxygen  reduction  at  the  cathode 

02  +  4e“  — 202"  (8) 


on  a  conventional  anode  (Ni-YSZ)  have  been  proposed  [13-17]. 
The  rate  expression  for  H2  oxidation  at  the  anode  used  in  this 
work  is  adopted  from  [15]  and  is  given  by 


RT 

^  p  ^tk ,  an 


f  2F?7H2,an 
V  RT 


—  exp 


^7H2,an 

RT 


(13) 


where  sn2,an  and  ^H2,an  are  the  anode  interface  conductivity 
(£2-1  m~2)  and  activation  overpotential  (V)  for  H2  oxidation, 
respectively.  su2m  is  estimated  using  the  following  empirical 
correlation: 


^H2,an  =  6.2  x  1011  x  exp 


— 120000  \  /  pu2o  \ 

RT  J  \  Req,H2PH2  J 


0.266 


(14) 


where  /Teq,H2  is  the  equilibrium  constants  for  the  reaction 
H2  +  I/2O2  — >  H2O.  Substitution  of  Eq.  (14)  into  Eq.  (13)  yields 
the  following  rate  expression  for  H2  oxidation: 


7H2,an  — 


2.1  x  1011  x  exp(-120000 /RT)p^RT 


(Keq,n2PH2f266F 


X 


exp 


(2Frm 2, 

V  RT 


an 


exp 


~Fm2, 


an 


RT 


(15) 


where  ?7H?,an  is  defined  as 

^?H2,an  —  0an  —  0e(^>  £a)  —  ^H2,an  (16) 


varies  depending  on  the  cathode  and  electrolyte  materials  as 
well  as  on  the  temperature.  For  conventional  cathode  material 
(LSM),  the  dissociation  of  adsorbed  oxygen  molecules  is  the 
rate-determining  step  (Minh  and  Takahashi,  1995).  Therefore, 
the  rate  expression  for  oxygen  reduction  at  an  SOFC  cathode  is 
based  on  this  step. 

The  rate  of  O2  reduction  at  the  cathode  is  determined  by  the 
rate  of  O2  adsorption  as  the  rate-determining  step  [14]  and  is 
given  by  [15] 


0.25  x  1010  x  exp(  —  1 30000 /RT)R 7/?Q^at 

F 


where  po2iCat  is  the  partial  pressure  of  oxygen  at  the  cathode. 
Ocat  is  the  activation  overpotential  for  O2  reduction  defined  as 


heat  —  0cat  0e(*3  £c)  ^cat 


(10) 


where  0Cat  and  0e(k  zc)  represent  the  cathode  and  electrolyte 
potentials,  respectively,  at  the  cathode/electrolyte  interface.  E®at 
is  the  reversible  potential  for  O2  reduction: 

kat  =  In  ^02, cat  (11) 


3.3.2.  H2  and  CO  oxidation  reactions 

Various  mathematical  expressions  for  the  rate  of  H2  oxidation 

H2  +  O2-  -*  H2O  +  2e"  (12) 


and  0an  and  (pQ{r ,  z&)  represent  the  anode  and  electrolyte  poten¬ 
tials,  respectively,  at  the  anode/electrolyte  interface.  an  is 
the  reversible  potential  for  H2  oxidation: 


E 


0 

H2,an 


TL)  in  (  ) 

2 FJ  \  ^eq,H2  PH2  J 


(17) 


Although  the  reaction  mechanisms  and  rate  expression  of  CO 
oxidation 


CO  +  O2-  -*  C02  +  2e“ 


(18) 


have  been  less  studied  than  that  of  H2  oxidation,  it  is  believed 
that  the  rate-determining  step  for  CO  oxidation  is  similar  to  that 
for  H2  oxidation  [13].  Setoguchi  et  al.  [18]  found  that  the  anodic 
interface  conductivity  for  CO  oxidation  (sco,an)  is  the  same  as 
that  of  H2  oxidation  if  the  partial  pressure  of  O2  during  H2 
oxidation  is  similar  to  that  during  CO  oxidation.  This  suggests 
that  the  interface  conductivity  of  CO  oxidation  can  be  described 
by  a  similar  correlation  to  that  of  H2  oxidation,  i.e., 


•^CO, an  =  6.2  X  1011  X  exp 


-1 20000  \  /  pc o2  \ 
RT  J  \  Keq  coPco  J 


0.266 


(19) 


where  7feq,co  is  the  equilibrium  constants  for  the  reaction 
CO  +  l/202  — >  C02. 

However,  it  was  recently  found  that  CO  oxidation  has  a  much 
higher  overpotential  than  H2  oxidation  for  the  same  current  den¬ 
sity  and  the  difference  cannot  be  attributed  to  differences  in 
concentration  overpotential  [19].  This  observation  agrees  with 
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that  of  Matzusaki  et  al.  [20]  who  found  the  rate  of  CO  oxida¬ 
tion  to  be  2-3  times  less  than  that  of  H2  oxidation  over  the  entire 
range  of  overpotential  and  the  same  partial  pressure  of  O2.  Based 
on  this  argument,  the  rate  of  CO  oxidation  reaction  used  in  this 
model  development  is  assumed  to  be  2.5  times  less  than  that  for 
H2  oxidation  reaction,  i.e.,  /co,an  =  0. 4  7h2,  an- Accordingly,  the 
rate  expression  for  CO  oxidation  is  expressed  as 


2cO,an  — 


0.84  x  1011  x  exp(  — 1 20000 /  RT) p/f/ R T 


X 


exp 


(K^coPCo)°-266F 
/  2Frjco  ,an  ___  f-Fri co  ,an 

v  RT 


exp 


RT 


(20) 


where  77co,an  is  defined  as 

0CO,an  =  0an  —  0e(G  Za)  —  ^CO,an  (21) 

and  an  is  the  ideal  anode  half-cell  potential  for  CO  oxidation 

reaction  as  given  by 


E 


o 

CO,  an 


f  Pco2  A 

\  ^eq,CO  PCO  / 


(22) 


3.4.  Chemical  reactions 

3.4.1.  Water-gas  shift  reaction  (WGSR) 

Due  to  its  high  operating  temperature,  nickel  can  act  as  a 
catalyst  for  the  WGSR  at  the  anode  when  syngas  is  used  as 
a  fuel  source.  The  rate  expression  Rwgsr  for  this  reaction  is 
written  as  [20,21] 

^wgsr  =  ^wgsr  (  Pu2o  Pco  -  ^H2^c°2  \  (23) 

V  Aeq,WGSR/ 


3.4.2.  Carbon  formation  reaction 

Carbon  formation  is  undesirable  since  it  can  impede  gas  flow 
and  block  active  sites  on  the  anode.  In  addition,  it  can  result  in  the 
growth  of  carbon  filaments  that  generate  massive  forces  within 
the  electrode  structure  which  can  cause  rapid  catalyst  break¬ 
down  [26].  Therefore,  it  is  important  to  understand  the  kinetics 
of  carbon  formation  and  determine  the  conditions  to  suppress 
its  formation.  For  syngas-fed  SOFCs,  carbon  formation  most 
likely  occurs  through  the  Boudouard  reaction,  as  indicated  in 
Eq.  (1),  since  Ni  is  a  good  catalyst  for  this  reaction.  The  reaction 
mechanism  and  kinetic  expression  for  the  Boudouard  reaction  on 
Ni-supported  YSZ  have  not  been  well  established.  However,  it 
is  believed  that  this  reaction  is  fast  enough  to  attain  equilibrium. 
Accordingly,  in  previous  studies,  a  thermodynamic  analysis  of 
carbon  formation  rather  than  a  kinetic  analysis  has  been  con¬ 
sidered  to  determine  suitable  operating  conditions  such  as  feed 
composition,  temperature  and  pressure  to  suppress  carbon  depo¬ 
sition  [14,26-29].  We  adopt  a  similar  approach  in  this  study  and 
do  not  incorporate  carbon  formation  directly  into  the  model  for 
cell  performance. 

Based  on  the  thermodynamics  of  the  Boudouard  reaction, 
carbon  can  form  when 

«C  =  ^eq.B/'co//'C02  >  1  (26) 

where  is  the  equilibrium  constant  of  the  Boudouard  reac¬ 
tion,  which  decreases  as  the  operating  temperature  rises.  The 
ratio,  ac ,  is  sometimes  referred  to,  more  or  less  correctly,  as 
“carbon  activity”.  The  possibility  of  carbon  formation  can  be 
obtained  from  the  solution  of  the  fuel  cell  model  that  determines 
values  of  pco  and  pco2  at  positions  along  the  anode.  The  value 
of  ac  at  each  position  is  then  assessed  according  to  the  crite¬ 
rion  given  in  Eq.  (26)  to  determine  whether  carbon  formation  is 
thermodynamically  possible. 


where  kwGSR  is  the  rate  constant  and  Ke q,wGSR  the  equilibrium 
constant  for  the  WGSR.  The  following  expression  for  kwGSR 
can  be  obtained  by  fitting  an  Arrhenius  curve  to  the  experimental 
data  given  in  [21]: 


o 

Mvgsr  =  1-21  x  10  x  exp 


103191 \ 
RT  ) 


(24) 


The  temperature  dependence  of  /%q,wGSR  is  given  by 

^eq,WGSR  =  exp(— 0.2935Z3  +  0.6351T3 

+  4.1788Z  +  0.3169)  (25) 

where  Z  =  2222  —  l 

Since  the  rate  of  the  WGSR  is  quite  high,  the  WGSR  has  been 
assumed  be  at  thermodynamic  equilibrium  in  various  SOFC 
models  [13,21-24].  However,  Ahmed  and  Foger  [25]  reported 
that  the  WGSR  reaches  very  close  to  equilibrium  only  at  high 
levels  of  fuel  utilization  (C/f  ~  80-100%).  The  actual  extent  of 
reaction  was  only  about  80-90%  of  the  equilibrium  level  at  lower 
fuel  utilization  (t/f  ~  0-0.80).  Therefore,  in  this  model,  we  do 
not  assume  that  the  WGSR  has  reached  equilibrium  and  instead 
consider  that  its  rate  is  given  by  Eq.  (23). 


3.5.  Boundary  conditions 

In  order  to  solve  the  system  of  equations,  boundary 
conditions  are  specified  at  all  the  outer  interfaces  of  the  com¬ 
putational  domains,  the  anode/electrolyte  (0  <  r<  rb,  z  =  za),  the 
cathode/electrolyte  interfaces  (0<r<r^,  z  =  zc)  and  at  the  anode 
tube  wall  (r=ra,  0<z<zt)-  At  the  fuel  inlet  (0<r<ra,  z  =  0), 
the  inlet  velocity  and  fuel  gas  mole  fractions  are  prescribed 
while  the  inlet  mole  fraction  of  O2  is  specified  at  the  air  inlet 
(0<r<rb,  z  =  L).  At  the  fuel  outlet  (ra<r<rb,  z  =  0),  only 
the  operating  pressure  is  prescribed,  i.e.,  pf=  1.013  x  105  Pa; 
for  all  other  variables,  the  gradients  in  the  flow  direction  are 
assumed  to  be  zero,  i.e.,  dvf/dz  =  0,  duf/dz  =  0,  dcfdz  =  0. 
Along  the  symmetry  line  at  r  =  0,  the  following  conditions 
hold:  dpf/dr  =  0 ,  dvf/dr  =  0,  duf/dr  =  0  and  dcfdr  =  0.  At  the 
anode/electrolyte  interface  at  z  =  za,  the  fuel  velocities  in  the 
radial  and  vertical  directions  are  zero,  i.e.,  uf  =  Vf  =  0.  The 
flux  of  each  gas  species  is  related  to  the  rate  of  H2  or  CO 

oxidation  according  to  Faraday’s  law  and  the  WGSR  kinet- 
— ►  — ► 

ics:  Ah2  =  -Dh23ch2/9z  =  Rwgsr  «/H2,an/2F,  Nh2o  = 
2H2,an/2 F  -  Rwgsr,  Nqo  =  -2co,an/2F  -  Rwgsr,  and 
Nqo2  =  2co,an/2F  +  Rwgsr-  The  charge  flux  is  related  to  the 
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(a)  Current  density  (A  cm-2)  (b)  Current  density  (A  cm-2) 

Fig.  4.  Experimental  (symbols)  and  predicted  (lines)  cell  performances  at  various  H2/N2  (a)  and  CO/N2  (b)  compositions  at  800  and  900  °C  for  a  fuel  flow  rate  of 
180  ml  min-1. 


(a)  Current  density  (A  cm'2)  (b)  Current  density  (A  cm'2) 


Fig.  5.  Experimental  (symbols)  and  predicted  (lines)  cell  performances  at  various  H2  +  CO2  (a)  and  CO  +  CO2  (b)  compositions  at  800  and  900  °C  for  a  fuel  flow 
rate  of  1 80  ml  min- 1 . 
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total  current  density  produced  by  H2  and  CO  oxidation,  i.e., 
— > 

i  =  —  cre30e/9^  =  Jh2  ,  an  +  7co,an-  At  the  cathode/electrolyte 
interface  at  z  =  zc,  the  molar  flux  of  O2  is  given  by 
Nq2  —  —Jca.t/4-F  and  the  current  density  i  =  —  /cat. 


3.6.  Model  calibration  and  numerical  implementation 


Before  the  cell  performance  could  be  modeled  and  validated 
against  experimental  data,  it  is  necessary  to  account  for  the  effect 
of  the  contact  resistance  at  the  current  collector/electrode  and 
electrode/electrolyte  interfaces.  The  following  approach  is  used. 
For  each  simulation,  the  cell  potential,  0C at  —  0an,  is  fixed  to 
an  arbitrary  value  and  the  model  equations  solved  to  yield  the 
current  density  distribution  along  each  electrode.  The  average 
current  density  along  each  electrode  is  then  determined  from 
this  distribution,  i.e., 


rrb 

/  ^(Al2,an  T  *^CO,an)dr 

JO 


(27) 


4.1.1.  Effect  of  fuel  flow  rate 

The  comparison  between  model  and  experimental  data  using 
pure  H2  as  fuel  at  various  flow  rates  for  cell  operation  at  800 
and  900  °C  is  presented  in  Fig.  3.  The  fitting  procedure  yields 
^contact  of  1.6  and  1.2  £2  cm2  at  800  and  900  °C,  respectively, 
which  are  in  the  range  reported  by  Jiang  [31]. 

As  evident  from  Fig.  3,  the  cell  potential  during  operation 
at  900  °C  is  slightly  greater  than  that  obtained  at  800  °C.  This 
occurs  because  the  open-circuit  is  higher  and  the  overpotentials 
are  lower  at  the  higher  temperature.  The  model  curves  slightly 
underestimate  the  experimental  data  at  800  °C  while  they  fit 
quite  well  at  900  °C.  They  also  indicate  that  cell  performance 
is  independent  of  fuel  flow  rate,  whereas  the  measured  values 
show  a  small  enhancement  in  cell  performance  at  higher  flow 
rates.  This  effect  is  more  apparent  at  the  lower  temperature.  For 
instance,  at  800  °C  and  0.05  A  cm-2,  the  cell  potential  increases 
from  ~0.86  to  ~0.9  V  when  the  fuel  flow  rate  is  raised  from  100 
to  180  ml  min-1 .  A  similar  effect  was  reported  by  Cunningham 
et  al.  [32]. 


Since  the  model  does  not  explicitly  account  for  the  contact  resis¬ 
tances  ^contact  in  the  cell,  the  actual  potential,  Vdelivered,  delivered 
to  the  external  load  is  estimated  as 

^delivered  =  (0  cat  —  0an)  —  7avg  ^contact  (28) 

The  value  of  ^contact  is  determined  by  calibrating  the  model  to 
experimental  data  using  the  following  procedure.  A  cell  perfor¬ 
mance  curve  at  a  given  set  of  gas  flow  rates  was  determined 
by  measuring  ^delivered  to  an  external  load  at  different  set  val¬ 
ues  of  the  load  resistance.  This  permitted  the  corresponding  7avg 
over  the  SOFC  to  be  obtained  and  a  performance  plot  of  Vdelivered 
versus  /avg  to  be  drawn.  The  model  was  calibrated  against  exper¬ 
imental  data  using  pure  H2  as  the  fuel  at  temperatures  of  800 
and  900  °C.  The  value  of  ^contact  was  obtained  by  fitting  the 
^delivered  —  4vg  plot  computed  according  to  Eq.  (28)  to  the  exper¬ 
imentally  observed  one  for  the  same  operating  conditions.  Since 
^contact  varies  only  with  temperature,  the  values  obtained  at  the 
two  temperatures  are  then  used  in  the  subsequent  simulations  at 
the  various  syngas  compositions. 

The  finite  element  software,  FEMLAB®  version  3.0  [30] 
was  employed  to  solve  the  system  of  equations.  The  num¬ 
ber  of  mesh  elements  in  solving  this  model  are  around  4000 
elements.  Depending  on  the  initial  guesses,  each  simulation 
required  between  2  and  5  min  for  completion  on  a  Pentium  IV 
1.4  GHz  with  1GB  RAM. 

4.  Results  and  discussion 

4.1.  Experimental  results  and  model  validation 

The  cell  performances  obtained  from  the  model  for  various 
syngas  compositions  are  compared  with  the  experimental  data 
at  800  and  900  °C  first  to  obtain  ^contact  and  then  to  validate  the 
model  for  other  operating  conditions.  The  model  input  parame¬ 
ters  and  operating  conditions  used  in  the  model  are  summarized 
in  Table  2.  The  cell  potential  was  varied  from  0.5  to  1.0  V  by 
setting  0cat  =  0  and  varying  0an  from  —0.5  to  —  1 .0  V. 


4.1.2.  Effect  of  N 2  dilution 

Fig.  4  shows  experimental  and  predicted  cell  performances 
with  H2  +  N2  and  CO  +  N2  mixtures  in  which  the  H2  (or  CO) 
composition  was  varied  from  ^100%  to  ~10%  at  both  800  and 
900  °C.  The  fuel  flow  rate  is  fixed  at  180  ml  min-1  in  all  cases. 
The  model  predictions  agree  well  with  the  experimental  data 
particularly  at  900  °C.  The  experimental  results  in  Fig.  4a  show 
that  the  cell  potential  improves  as  H2  concentration  increases 
although  the  OCPs  are  similar.  The  model  simulations  indicate 
that  one  reason  for  this  effect  is  that  N2  acts  as  an  inert  gas  to 
the  diffusion  of  H2  from  the  bulk  to  the  anode.  Consequently, 
H2  is  less  evenly  distributed  over  the  electrode.  In  addition, 
from  a  thermodynamic  point  of  view,  a  higher  partial  pressure 
of  H2  leads  to  a  higher  cell  potential.  However,  the  impact  of 
N2  dilution  is  modest  and  the  cell  performance  decreases  by 
about  15%  when  the  H2  concentration  is  reduced  from  ^100% 
to  ~10%. 

When  CO  is  diluted  with  N2  at  800  °C  (Fig.  4b),  the  OCP 
is  similar  to  that  for  the  H2  +  N2  mixture  with  the  same  fuel 
gas  composition,  but  the  cell  performance  is  lower.  On  the  other 
hand,  the  cell  performances  using  CO  +  N2  and  H2  +  N2  mixtures 
are  much  closer  to  each  other  at  900  °C. 

Considerable  discrepancies  between  the  model  results  and 
experimental  data  are  observed  at  800  °C  (Fig.  4b),  particularly 
at  high  current  density  and  low  CO  concentration.  The  overes¬ 
timation  of  cell  performance  by  the  model  is  likely  due  to  the 
effect  of  carbon  formation  via  the  Boudouard  reaction  that  is 
not  considered  in  the  model.  Carbon  formation  is,  of  course, 
possible  when  CO  is  used  as  fuel  and  most  likely  at  low  cur¬ 
rent  density  when  pco  remains  relatively  high  and  pco2  is  low. 
Each  of  the  experimental  performance  curves  in  this  study  was 
obtained  by  varying  the  current  density  from  low  to  high  values 
using  the  same  electrodes  without  disassembling  the  apparatus. 
Since  carbon  is  not  easily  removed  once  formed  [33],  it  is  likely 
that  some  carbon  formed  at  low  current  density  remained  on  the 
anode  when  higher  current  densities  were  applied  and  continued 
to  reduce  cell  performance. 
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Table  3 

Fuel  compositions  due  to  the  WGSR  at  900  °C  for  various  H2/CO2  mixtures 


Inlet  composition  (%) 

Fuel  composition  at  SOFC  anode  (%) 

E° 

cell 

h2 

C02 

h2 

CO 

C02 

h2o 

99 

0 

99 

0 

0 

1 

1.15 

80 

20 

63.3 

16.7 

3.3 

17.7 

0.98 

60 

40 

34.4 

25.6 

14.4 

26.6 

0.93 

40 

60 

14.4 

24.5 

34.4 

26.6 

0.89 

20 

80 

3.3 

16.7 

63.3 

17.7 

0.83 

4.1.3.  Effect  of  CO 2  dilution 

Fig.  5  shows  experimental  and  predicted  cell  performances 
using  H2  +  CO2  and  CO  +  CO2  mixtures  at  both  800  and  900  °C. 
The  composition  of  H2  and  CO  was  varied  over  the  range  from 
100%  to  10%  in  each  of  these  mixtures.  The  effect  of  CO2 
dilution  is  more  pronounced  than  that  of  N2  dilution.  The  cell 
performance  with  H2  (or  CO)  +  C02  is  lower  than  that  with 
H2  (or  CO)  +  N2  for  the  same  fuel  concentration.  The  OCP 
decreases  significantly  as  CO2  concentration  is  increased  when 
either  H2  or  CO  is  the  fuel.  As  shown  in  Fig.  5a,  the  OCP  at  800 
and  900  °C  decreases  by  20%  and  26%,  respectively  as  H2  con¬ 
centration  is  decreased  from  100%  to  20%.  This  is  likely  due  to 
the  impact  of  the  WGSR.  For  the  H2  +  CO2  mixture,  the  WGSR 
strongly  influences  the  equilibrium  compositions  by  converting 
some  of  H2  and  CO2  to  H2O  and  CO  which  leads  to  a  lower  OCP, 
as  indicated  by  Eqs.  (17)  and  (22).  The  effect  of  CO2  dilution 
on  the  OCP  at  900  °C  is  summarized  in  Table  3. 

From  Fig.  5b,  it  is  evident  that  at  high  CO  concentration 
(^40%  to  ~100%)  the  same  trends  of  experimental  and  model 
predicted  cell  performances  are  observed  at  both  800  and  900  °C. 
Unfortunately,  the  model  slightly  underestimates  the  experimen¬ 
tal  cell  performance  at  both  temperatures.  On  the  other  hand,  the 
experimental  cell  potentials  at  low  CO  concentrations  (~5%  to 
^20%)  show  a  drop  at  high  current  densities  that  is  not  predicted 
by  the  model  especially  at  800  °C.  These  discrepancies  are  also 
likely  caused  by  carbon  formation  on  the  anode. 

4.1.4.  Effect  ofH2/CO  composition 

Fig.  6  shows  the  experimental  and  predicted  cell  perfor¬ 
mances  for  H2/CO  ratios  diluted  with  70%  N2  to  reduce  the  rate 
of  carbon  formation  that  depends  on  the  CO  concentration  in  the 
fuel  gas.  The  model  only  slightly  underestimates  the  measured 
performance  data  at  both  temperatures.  The  OCP  is  essentially 
independent  of  the  relative  proportions  of  FU  and  CO.  Moreover, 
the  cell  performance  achieved  when  CO  is  the  only  oxidizable 
gas  is  slightly  lower  than  that  obtained  when  FU  is  the  only  fuel. 
This  means  that  even  though  the  CO  oxidation  rate  is  less  than 
the  H2  oxidation  rate,  it  has  only  a  small  impact  on  the  relative 
cell  performance.  The  main  reason  for  this  is  that  the  cathodic 
overpotential  plays  a  larger  role  on  the  cell  performance  than  the 
anodic  overpotential  [33]. 

4.1.5.  Effect  of  simulated  syngas  composition 

In  order  to  evaluate  the  ability  of  the  model  to  predict  cell 
performance  with  syngas  fuel,  model  predictions  were  compared 
to  experimental  results  obtained  using  various  gas  compositions. 
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Fig.  6.  Experimental  (symbols)  and  predicted  (lines)  cell  performances  at  vari¬ 
ous  H2/CO  compositions  diluted  with  70%  N2  at  (a)  800  °C  and  (b)  900  °C  for 
a  fuel  flow  rate  of  1 80  ml  min- 1 . 

The  gas  compositions  used  in  these  experiments  are  listed  in 
Table  4. 

Experimental  and  predicted  cell  performances  for  various 
syngas  compositions  at  800  and  900  °C  are  presented  in  Fig.  7. 
The  predictions  agree  well  with  the  experimental  data  at  900  °C, 
but  underestimate  the  data  at  800  °C  although  there  is  close 
agreement  with  regard  to  the  OCR 

Fuel  FI  consisting  of  moist  H2  provides  the  highest  perfor¬ 
mance.  However,  it  is  interesting  to  note  that  fuel  F6  which  is  a 
mixture  of  20%  H2,  20%  CO  and  60%  N2  performs  slightly  bet- 


Table  4 


Simulated  syngas  composition 


Fuel  no. 

Gas  composition 

h2 

h2o 

CO 

C02 

n2 

FI 

97 

3 

— 

— 

— 

F2 

20 

3 

— 

— 

80 

F3 

20 

3 

— 

14 

66 

F4 

20 

3 

20 

14 

43 

F5 

32 

3 

45 

15 

3 

F6 

20 

3 

20 

0 

57 

R.  Suwanwarangkul  et  al.  /  Journal  of  Power  Sources  161  (2006)  308-322 


317 


Fig.  7.  Experimental  (symbols)  and  predicted  (lines)  cell  performances  for  sim¬ 
ulated  syngases  at  compositions  summarized  in  Table  4  at  (a)  800  °C  and  (b) 
900  °  C  for  a  fuel  flow  rate  of  1 80  ml  min- 1 . 

ter  than  moist  H2  diluted  by  80%  N2  (F2).  This  indicates  that  CO 
should  be  a  useful  fuel  for  SOFCs.  Fuel  F3  is  a  mixture  of  H2, 
CO2  and  N2,  whereas  F4  and  F5  are  examples  of  syngas  com¬ 
positions  typically  obtained  from  air  and  oxygen  gasifiers  [34], 
respectively.  The  OCPs  obtained  using  F3,  F4  and  F5  are  less 
than  those  of  FI ,  F2  and  F6,  likely  due  to  the  effect  of  the  WGSR 
in  the  presence  of  CO2,  as  previously  described.  The  cell  per¬ 
formance  using  the  F5  fuel  is  slightly  greater  than  that  using  F4 
because  of  the  higher  fuel  content  and  lower  N2  dilution.  How¬ 
ever,  when  the  cost  of  syngas  production  is  considered,  the  use 


of  syngas  composition  F4  may  be  the  most  economical  choice. 
F3  yields  the  poorest  cell  performance  because  of  its  low  fuel 
content  and  considerable  dilution  by  presence  of  CO2  and  N2. 
If  CO2  is  removed  before  the  syngas  enters  the  SOFC,  the  cell 
performance  could  be  further  increased.  On  the  other  hand,  the 
presence  of  CO2  helps  prevent  carbon  formation  on  the  anode. 
Thus,  a  trade-off  must  be  made  to  maximize  cell  performance 
and  anode  durability. 

4.1.6.  Cell  inspection  after  operation 

Once  an  experiment  was  completed  at  each  operating  tem¬ 
perature,  the  test  rig  was  disassembled  and  the  electrochemical 
cell  visually  inspected.  A  significant  amount  of  Ni  and  carbon 
was  found  to  have  deposited  on  the  anode  current  collector  after 
operation  at  800  °C,  as  confirmed  by  XRD  analysis.  Decoloura¬ 
tion  of  the  anode  surface  indicating  Ni  depletion  was  also  noted. 
A  thick  layer  of  carbon  was  also  observed  in  the  anode  feed 
entrance  region  of  the  cell. 

On  the  other  hand,  only  a  small  part  of  the  Ni  anode  material 
delaminated  and  deposited  on  the  anode  current  collector  after 
cell  operation  at  900  °C.  No  evidence  of  carbon  was  found  at  the 
anode  current  collector  or  the  anode  surface  and  only  a  small 
amount  was  observed  at  the  anode  feed  entrance  region. 

SEM  images  of  anode  surfaces  obtained  after  operation  at  800 
and  900  °C  are  presented  in  Fig.  8.  The  anode  surface  obtained 
at  800  °C  appears  to  be  worse  than  that  obtained  at  900  °C.  The 
anode  surface  at  the  lower  temperature  has  become  almost  com¬ 
pletely  deactivated  due  to  the  presence  of  many  particles  that 
block  the  pores  and  prevent  reactant  gases  from  reaching  the 
active  sites  within  the  pores.  Those  particles  are  should  be  com¬ 
posed  of  Ni  delaminated  from  the  anode  surface.  Furthermore, 
the  pore  size  after  operation  at  800  °C  appears  smaller  than  that 
at  900  °C  likely  due  to  Ni  agglomeration. 

From  visual  and  SEM  investigation,  anode  delamination 
appears  to  be  related  to  carbon  formation  which  deteriorates 
cell  performance.  The  actual  mechanism  of  anode  delamination 
due  to  carbon  formation  has  not  been  clearly  established.  How¬ 
ever,  the  role  of  Ni  used  in  an  SOFC  anode  is  likely  similar  to 
that  in  Ni-supported  catalysts.  A  mechanism  for  Ni  disintegra¬ 
tion  from  its  support  due  to  carbon  formation  on  Ni-supported 
catalysts  has  been  proposed  by  Snoeck  et  al.  [35].  According  to 
this  explanation,  the  Boudouard  reaction  leads  to  the  adsorption 
of  carbon  atoms  onto  the  catalyst  which  then  dissolve  into  the 


Fig.  8.  SEM  images  of  anode  surface  obtained  after  operation  at  800  °C  (a)  and  900  °C  (b). 
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Fig.  9.  Surface  and  arrow  plots  showing  the  velocity  field  inside  fuel  channel  for  the  base  case  conditions. 


nickel  on  the  gas  side  and  diffuse  through  the  Ni  particle  to  the 
support  side,  which  eventually  causes  the  separation  of  Ni  from 
the  support  leading  to  nucleation  of  carbon  filaments. 

4.2.  Simulation  results  for  simulated  syngas 

The  simulation  results  discussed  here  are  based  on  simulated 
syngas  obtained  from  biomass  gasification  in  air  (fuel  F4  in 
previous  section).  The  operating  conditions  are  summarized  in 
Table  5.  It  is  noted  here  that  the  contact  resistance  is  assumed  to 
be  zero  for  the  purposes  of  this  simulation. 

The  velocity  field  inside  the  fuel  channel  for  the  base  case 
conditions  is  depicted  in  Fig.  9.  The  velocity  is  high  at  the 
entrance  region  and  is  reduced  as  one  moves  away  from  the 
entrance.  Since  the  feed  tube  is  much  smaller  than  the  fuel  chan¬ 
nel,  uneven  flow  occurs  as  the  fuel  gas  is  dispersed  from  the  feed 
tube  to  the  fuel  channel.  Once  the  fuel  gas  reaches  the  anode 
surface,  it  begins  to  flow  in  the  reverse  direction  toward  the  dis¬ 
charge  from  the  fuel  channel.  It  is  found  that  the  gas  velocity 
is  almost  zero  near  the  edge  of  the  fuel  channel  (rb  =  0.018  m). 
Accordingly,  transport  of  the  gases  in  this  region  occurs  primar¬ 
ily  by  diffusion  and  eventually  local  current  density  in  this  area 
as  illustrated  in  Fig.  10.  Therefore,  reactant  concentrations  (H2 
and  CO)  at  the  edge  of  the  fuel  channel  are  less  than  those  in  the 
fuel  channel  close  to  the  feed  tube. 

The  distributions  of  the  H2,  CO,  CO2,  H2O  concentrations 
inside  the  fuel  channel  are  displayed  in  Fig.  10.  The  H2  concen¬ 
tration  declines  sharply  as  the  gas  moves  through  the  fuel  inlet 
region  and  then  decreases  more  gradually  in  the  region  further 
away  from  the  fuel  inlet.  This  is  due  to  its  high  rate  of  consump- 

Table  5 


Base  case  operating  conditions 


Operating  condition 

Value 

Pressure 

1  atm  (101.35  xlO5  Pa) 

Temperature 

900  °C 

Fuel  flow  rate 

180  ml  min-1 

Fuel  gas  composition  (fuel,  F4) 

20%  H2,  20%  CO,  14%  C02 
and  3%  H20,  43%  N2 

Cathodic  half-cell  potential  (0C at) 

OV 

Anodic  half-cell  potential  (0an) 

-0.7  V 

Oxidant 

Air 

Contact  resistance  (/^contact) 

0  Q  cm-2 

tion  by  electrochemical  oxidation  and  the  WGSR  to  produce 
CO  and  H2O  on  the  portion  of  the  anode  near  the  inlet.  This  is 
confirmed  in  Fig.  10c  that  shows  the  CO  concentration  rising  at 
the  anode  surface  near  the  fuel  inlet.  Once  the  WGSR  reaches 
equilibrium,  H2  is  mainly  consumed  by  electrochemical  oxida¬ 
tion  which  is  slower  than  the  WGSR,  leading  to  a  more  uniform 
concentration.  The  H2O  concentration  increases  abruptly  at  the 
fuel  inlet  region  due  to  the  WGSR.  This  also  causes  the  concen¬ 
tration  profiles  of  CO  and  CO2  in  fuel  inlet  region  to  be  inverted 
with  respect  to  one  another. 

As  shown  in  Fig.  11,  the  oxygen  concentration  distribution 
remains  uniform  throughout  the  air  channel.  This  implies  that 
sufficient  amount  of  O2  is  being  transferred  to  the  cathode  sur¬ 
face  even  if  its  source  is  ambient  air.  Therefore,  the  use  of  an  air 
feed  tube  is  not  necessary  in  this  set-up. 

The  distributions  of  the  total  current  density  and  the  partial 
current  densities  due  to  H2  and  CO  oxidation  along  the  electrode 
surface  are  shown  in  Fig.  12.  These  show  that  the  current  density 
gradually  decreases  along  the  length  of  the  anode  surface.  The 
rate  of  H2  oxidation  remains  about  2.5  times  greater  than  that 
of  CO  oxidations,  as  specified  by  the  electrochemical  kinetic 
expressions,  although  the  diffusivity  of  CO  is  much  less  than 
that  of  H2  (Table  2).  This  indicates  that  the  influence  of  mass 
transport  along  the  anode  surface  is  not  significant. 

The  breakdown  of  the  overpotentials  along  both  electrodes  is 
shown  in  Fig.  13.  All  overpotentials  are  quite  uniform  along  the 
electrode  surface  due  to  the  uniform  current  density  distribution. 
The  cathodic  overpotential  (?7cat)  is  about  two  times  larger  than 
the  other  overpotentials.  The  ohmic  overpotential  (r/^),  anodic 
overpotential  for  H2  oxidation  (/7H?,an)  and  anodic  overpotential 
for  CO  oxidation  (r/cctan)  are  very  similar  in  magnitude.  Since 
the  cathodic  overpotential  is  large  and  not  dependent  on  the 
H2/CO  ratio,  the  cell  potential  is  not  affected  by  the  H2/CO 
ratio,  as  observed  in  Fig.  6. 

4.3.  Thermodynamic  analysis  of  carbon  formation  on  an 
SOFC  anode 

As  discussed  previously,  the  possibility  of  carbon  for¬ 
mation  is  based  on  whether  the  carbon  activity  ratio  ac  = 
^eq,B  Pco/ Pco2  obtained  from  the  model  exceeds  1 .0  (Eq.  (26)). 
This  calculation  accounts  for  the  influences  of  the  electrochem¬ 
ical  oxidation  of  H2  and  CO,  operating  temperature,  pressure 


R.  Suwanwarangkul  et  al.  /  Journal  of  Power  Sources  161  (2006)  308-322 


319 


Max:  0.2 


0.005 

0.004 

gT  0.003 
N  0.002 

0.001 

0.0 


n 


0.19 
0.18 
0.17 
0.16 

Min:  0.151 


Fuel  outlet 

i - 1 - 1 - 1 — i ; r 

0.003  0.006  0.009 


•  i  • 


1  1  1  *  i  1  '  1  T 

0.012  0.015  0.018 


r  (m) 


0.003 


0.006 


0.009 
r  (m) 


0.012 


0.015 


Min:  0.20 


0.018 

Max:  0.14 


0.005 

0.004 

9  0.003 
SI  0.002 

0.001 

0.0 


Fuel  inlet 
t - 1 - : - 


Fuel  outlet 


0.0 


0.003 


1 — I-7- 

0.006 


I  i  i 


0.009 

r(m) 


0.012 


0.015 


0.018 


Fig.  10.  Surface  and  contour  plots  of  concentrations  for  H2  (a);  H2O  (b);  CO  (c)  and  CO2  (d)  inside  fuel  channel. 


and  inlet  syngas  composition.  ^contact  has  been  set  to  zero  for 
the  purposes  of  these  calculations.  The  effect  of  the  CO  content 
and  current  density  on  the  carbon  activity  ratio  for  cell  operation 
at  800  and  900  °C  is  depicted  in  Fig.  14.  The  value  of  the  carbon 
activity  ratio  shown  is  the  average  of  the  local  values  computed 
along  the  anode  surface,  i.e., 


^c,ang 


‘>'b 


rac(r)  dr 


(29) 


The  CO  content  is  defined  as 

mole  fraction  of  CO 

CO  content  =  - 

mole  fraction  of  CO  +  mole  fraction  of  H2 

(30) 

The  inlet  fuel  stream  is  assumed  to  contain  no  H2O  and  CO2. 

Fig.  14  indicates  that  the  CO  content  and  current  density  have 
a  significant  impact  on  the  possibility  of  carbon  formation  (i.e., 
ac  >  1)  for  cell  operation  at  800  °C.  In  contrast,  carbon  forma- 
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Fig.  11.  Surface  and  contour  plots  of  O2  concentration  inside  air  channel. 


Fig.  12.  Current  density  distributions  along  electrode  surface. 


tion  can  be  avoided  at  900  °C  over  the  entire  range  of  current 
density  and  CO  contents  because  of  the  lower  equilibrium  con¬ 
stant  for  the  Boudouard  reaction  at  the  higher  temperature.  This 
difference  agrees  with  the  experimental  observations  described 
in  Section  4.1.6.  Nonetheless,  operation  at  800  °C  or  below  is 
potentially  more  attractive  because  of  the  wider  choice  of  mate- 


Fig.  14.  Effects  of  current  density  and  CO  content  on  carbon  activity  at  1  atm 
and  800  and  900  °C. 


Fig.  13.  Break  down  of  overpotentials  along  electrode  surface. 


Fig.  15.  Effects  of  current  density  and  CO  content  on  carbon  activity  at  900  °C 
and  various  operating  pressures. 
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Fig.  16.  Effects  of  current  density  and  CO  content  on  carbon  activity  at  800  °C, 
3  atm  and  various  H2O  compositions. 

rials  that  can  be  used  for  cell  components.  It  is  also  observed 
that  likelihood  of  carbon  formation  at  800  °C  can  be  reduced 
by  operating  at  higher  current  density  and  lower  CO  content. 
An  increase  in  current  density  enhances  the  rate  of  H2  and  CO 
oxidation  and  the  rate  of  H2O  and  CO2  production.  Thereafter, 
H2O  reacts  with  CO  via  the  WGSR  to  form  more  CO2  and  H2. 
These  effects  tend  to  reduce  at  lower  ac- 

In  actual  SOFC  operation,  syngas  fuel  typically  contains 
about  50-60%  CO.  Therefore,  according  to  Fig.  14,  the  cell 
should  be  operated  at  a  minimum  of  0.1  A  cm-2  to  ensure  that 
enough  CO2  and  H2O  are  generated  to  prevent  carbon  formation 
at  800  °C.  However,  it  should  be  noted  that  this  calculation  is 
based  on  a  syngas  containing  only  H2  and  CO.  If  the  inlet  syn¬ 
gas  composition  contains  H2O  and/or  CO2,  operation  of  a  lower 
current  density  than  0. 1  A  cm-2  might  be  possible  because  their 
presence  helps  to  prevent  carbon  formation,  as  will  be  described 
later. 

Fig.  15  displays  the  effects  of  current  density  and  CO  content 
on  the  carbon  activity  ratio  at  various  operating  pressures  and 
900  °C.  An  increase  in  operating  pressure  favours  carbon  for¬ 
mation,  especially  at  high  CO  content  and  low  current  density. 
Operation  at  3  atm  has  the  danger  of  forming  carbon  over  the 
entire  range  of  operating  current  densities  and  CO  contents.  The 
application  of  such  a  pressure  has  recently  received  attention 
due  to  the  proposal  for  the  use  a  SOFC/micro-turbine  hybrid 
cycle  which  has  a  high  overall  efficiency,  but  requires  a  high 
inlet  pressure  for  the  micro-turbine. 

As  previously  described,  an  operating  target  of  800  °C  and 
3  atm  would  be  ideal.  Unfortunately,  these  conditions  can  lead 
to  unreliable  and  low  cell  performance  due  to  carbon  forma¬ 
tion.  One  of  the  possible  ways  to  reduce  this  likelihood  is  to 
increase  the  amount  of  H2O  and  CO2  in  the  inlet  stream.  The 
effects  of  the  current  density  and  CO  content  on  a  at  various 
inlet  H2O  and  CO2  compositions  at  800  °C  and  3  atm  are  shown 
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Fig.  17.  Effects  of  current  density  and  CO  content  on  carbon  activity  at  800  °C 
and  3  atm  and  various  CO2  compositions. 

in  Figs.  16  and  17,  respectively.  The  likelihood  of  carbon  for¬ 
mation  is  substantially  decreased  by  increasing  H2O  or  CO2 
concentration  in  the  inlet  stream  especially  at  high  CO  content. 
A  syngas  fuel  containing  50-60%  CO  content  requires  at  least 
20%  H2O  or  CO2  concentration  to  reduce  the  risk  of  carbon  for¬ 
mation  during  operation.  Unfortunately,  cell  power  will  slightly 
be  reduced  since  this  tends  to  lower  the  OCP. 

5.  Conclusions 

A  mechanistic  model  based  on  a  button  cell  geometry  is  vali¬ 
dated  against  experimental  performance  data  at  800  and  900  °C. 
This  comparison  shows  that  the  model  agrees  reasonably  well 
with  the  experimental  data  at  900  °C  but  not  as  well  in  the  case 
of  800 °C. 

The  experimental  data  and  the  model  simulations  indicate 
that  the  major  reason  for  achieving  lower  cell  performance  when 
the  SOFC  operates  with  syngas  fuel  is  carbon  formation  via 
the  Boudouard  reaction  that  damages  the  anode  surface.  If  this 
problem  can  be  alleviated,  then  the  analysis  shows  that  CO  is 
an  effective  fuel.  CO2  dilution  has  a  more  pronounced  effect  on 
the  cell  performance  than  that  of  N2  because  CO2  significantly 
reduces  OCP.  A  thermodynamic  analysis  of  carbon  formation 
shows  that  the  operating  pressure,  temperature  and  inlet  H2  and 
CO  compositions  are  important  factors.  The  operation  of  SOFCs 
at  intermediate  temperature  and  high  pressure  with  syngas  is 
limited  by  carbon  formation.  However,  this  may  be  overcome 
by  adding  H2O  or  CO2  into  the  syngas  fuel,  although  at  the 
expense  of  some  loss  in  the  electric  power  delivered  by  the  fuel 
cell. 
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